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We establish new Bonferroni-type lower bounds for the probability of a 
union of finitely many events where the selection of intersections in the 
estimates is determined by the clique complex of a chordal graph. 

1 Introduction 



The classical Bonferroni inequalities state that for any finite collection of events {A 1 ,} u6 y, 
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where £P*(V) denotes the set of non-empty subsets of V. Numerous variants of these 
inequalities are known, see, e.g., Galambos and Simonelli [6] for a survey. 

In this paper, we establish a new variant where the selection of intersections in the 
probability bounds is determined by a chordal graph. Recall that a graph is chordal if 
it contains no cycles of length four or any higher length as induced subgraphs. 

We refer to [2] for terminology on graphs. We write G = (V,E) to denote that G is a 
graph having vertex-set V, which we assume to be finite, and edge-set E, consisting of 
two-element subsets of V. We use ^(G) to denote the clique complex of G, that is, the 
abstract simplicial complex consisting of all non-empty cliques of G. 

A recent connection between Bonferroni inequalities and chordal graphs is the chor- 
dal graph sieve, which states that in the classical Bonferroni upper bound, £P*(V) may 
be replaced by the clique complex of some arbitrary chordal graph on V. We recall the 
chordal graph sieve in the following proposition. 
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Proposition 1.1 ([3j[4j). Let {A v } veV be a finite collection of events, where the indices form 
the vertices of a chordal graph G. Then, 



and 



Pr 



\JA V < £ (-l)HMpr f|A 

\veV ) ietf(G) 



(2) 



iel ) 



Pr 



J A, < £ (-l)I^Pr Qa, 



fe^G) 
|7|<2r-l 



(r = 1,2,3,...). 



(3) 



V iei / 



Remark 1.2. Inequality ([2) interpolates between Boole's inequality (G edgeless) and the 
sieve formula (G complete). Moreover, it generalizes a well-known result of Hunter 
IH on trees. Note that for any tree G = (V,E), the clique complex of G consists of the 
vertices and edges of G. 

Remark 1.3. The case r = 2 in (3j was independently found by Boros and Veneziani [T| 
by linear programming techniques; see also [1161] for a discussion on this bound. 



2 Main result 



In the sequel, we use c(G) resp. a(G) to denote the number of connected components 
of G resp. the independence number of G. Our main result, which is a lower bound 
analogue of Proposition [TTTJ reads as follows. 

Theorem 2.1. Let {A v } veV be a non-empty finite collection of events, where the indices form 
the vertices of a chordal graph G. Then, 
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Remark 2.2. As an immediate consequence of (O, and in contrast to the classical Bon- 
ferroni inequalities, the lower bound in © is guaranteed to be non-negative. 

Remark 2.3. Inequality © becomes an equality if G is complete. In this case, |4| agrees 
with the sieve formula, while (O coincides with the usual Bonferroni lower bounds. 

The proof of Theorem 12.11 is facilitated by a lemma and a proposition. The lemma 
contains a well-known combinatorial identity. 
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Lemma 2.4 ([13]]). For any n e N and any m e N , 

Proposition 2.5. For any chordal graph G = ( V ,E), 

(-if 1 ' 1 < c(G) (r = 1,2,3,...), (7) 

|J|<2r 

with equality if r > \V\/2. 

Proof. It suffices to show that {7J holds for any connected chordal graph G, that is, 

J^(-lf-i<l ( r = 1,2,3,...), (8) 

|J|<2r 

with equality if r > |V|/2. If G has at most one vertex, then this statement trivially 
holds. We proceed by induction on the number of vertices of G. If G has at least two 
vertices, then since G is chordal, G contains a simplicial vertex s, that is, a vertex whose 
neighborhood Ng(s) is a clique. Since G is connected, Ng(s) * 0. Moreover, since s is 
simplicial, the subgraph G - s obtained from G by removing s is connected. Since, in 
addition, G - s is chordal, we may apply the induction hypothesis to G - s, whence 
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where the equality in ([9]) follows from {6|. Note that if r > \ V\/2, then by the induction 
hypothesis, the first inequality in ([9j becomes an identity. The same applies to the 
second inequality, since for r > \V\/2 the binomial coefficient ( 2r-i ) vanishes. □ 

Remark 2.6. Eq. (HJ can also be proved by combining the following topological results: 

1. The clique complex of any connected chordal graph is contractible (|E|). 

2. For any contractible abstract simplicial complex 5? ', the Euler characteristic of its 
(2r- l)-skeleton is at most 1, with equality if r > \V\/2 (||14j|). 



We are now ready to proceed with the proof of Theorem \2A 
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Proof of Theorem \2.1\ It suffices to prove since ([4| follows from |[5| by choosing some 
arbitrary r > \ V\/2. For any I C V, I * 0, define 



(10) 



Note that the Bj's form a partition of {J ve yA v , that is, 

\Ja v = |J B, and 1^0.8, = (J*/). 



Clearly G[/] is chordal and c(G[/]) < a(G) for any / C V, J * 0. Hence, by applying 
Proposition [23] to G[/] we obtain 
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Remark 2.7. In view of the preceding proof, the bounds can be sharpened by replacing 
a(G) with a'(G) := max{c(G[/]) | / c V, J * 0, Bj * 0} where Bj is defined as in {T0|. In 
particular, if the underlying graph matches the structure of the events in such a way 
that G[J] is connected whenever Bj is non-empty, then a'(G) = 1 and the bounds in {2| 
and (HI coincide. The resulting identity is well-known in abstract tube theory |3|. 

Remark 2.8. The requirement that G is chordal cannot be omitted from Theorem 12. II 
Consider the non-chordal graph G = (V,E) depicted in Figured] and consider events 
A v , v € V, with Pi(A v ) = 1 for any v € V. The clique complex of this graph consists of 8 
cliques of size 1, 20 cliques of size 2, and 16 cliques of size 3. Since a(G) = 3, the first 
inequality in Theorem 12 . 1 1 gives the non- valid bound 1 > ±(8 - 20 + 16) = f . 

Remark 2.9. The graph G in the preceding counterexample was obtained as a subgraph 
of G 3 — the first graph in an infinite sequence (Gfc)fc=3,5,... of non-chordal graphs for 
which the hypothesis of Theorem 12.11 does not hold. Each G^ in this sequence is the 
join of k disjoint copies of the edgeless graph on three vertices. It turns out that G^ 
has independence number 3, and that the clique complex of G/ c has Euler characteristic 
1 + 2 k . By considering events A v , v e V, with Pr(A„) = 1 for any vertex v of G^, the first 
inequality in Theorem l2. 1 I becomes 1 > ^(1 + 2 k ), which is false for k = 3, 5, 
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Figure 1: A non-chordal graph for which the hypothesis of Theorem [2J] does not hold. 



Theorem 2.10. Under the requirements of Theorem \2.1\ 
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Proof. In view of the proof of Theorem l2.1[ it suffices to show that for any / C V, } 0, 
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A two-fold application of Proposition 12.51 reveals that the sum on the left-hand side 
equals c(G[/]), while the sum on the right-hand side is bounded above by c(G[/]). □ 



3 Particular cases 

In case of a tree, Theorem 12 . 1 1 gives a lower bound analogue of Hunter's inequality (9|- 
Hunter's upper bound on Pr(\J veV A v ) coincides with the bracketed term in (fTTT) . 

Corollary 3.1. Let {A v } veV be a non-empty finite collection of events, where the indices form 
the vertices of a tree G = (V,E). Then, 



Pr 



[jA v 
\ veV J 



> 



a(G) 



jyr(A v )- Pr(A v nA, 

\ veV [v,w}eE 



(11) 



Remark 3.2. Since a(G) < n - 1 for any tree G on n > 1 vertices, we conclude from (II 1 1> 
that, in this case, 
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To obtain a maximum lower bound of this type, we adapt Hunter's method |9| for 
computing a minimum upper bound: Any minimum spanning tree of the complete 
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weighted graph on V with edge weights Pr(A v C\A W ) for any v,w € V, v * w, maximizes 
the right-hand side of the preceding inequality. Such a minimum spanning tree can be 
found in polynomial time, e.g., by applying Kruskal's algorithm [I10IJ . 

Remark 3.3. If G is a path on n vertices, then a(G) = \n/2] and hence, 
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In order to find a path on V that maximizes the lower bound in (|12H , consider the com- 
plete graph on V with edge lengths Pr(Aj, (~)A W ) for any v,u> € V, v * w. Any minimum 
length Hamiltonian path in this graph maximizes the right-hand side of $12) . Unfortu- 
nately the problem of finding a minimum length Hamiltonian path is NP-hard [8|. 



Corollary 3.4. Let A\,...,A n be anon-empty finite collection of events. Then, 
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Proof. The inequality follows from 1)12)) by averaging over all paths on V = {1, . . ., n). □ 

Remark 3.5. Inequality (113ft may be considered as a lower bound analogue of KwerePs 
inequality HU, which states that Pr((J ve y A v ) < £? =1 Pr(A,) - § ^ 1 < l<y <„ Pr(A, n A ; ). 

Our next corollary provides a lower bound analogue of Seneta's inequality [I15IJ . The 
upper bound in Seneta's inequality coincides with the bracketed term that follows. 

Corollary 3.6. Let Ai,...,A n be a finite collection of events, and j,k € {!,..., n}. Then, 
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provided n> 2 - b^ where 5^ denotes the Kronecker delta. 



Proof. Define G = K 2 s . k *L n _ 2+ t, jk (the join of a complete graph on 2 - b^ vertices and 
an edgeless graph on n - 2 + b& vertices), and apply Theorem l2.ll □ 

By taking the average over all distinct j,k = l,...,n in the preceding inequality, a 
lower bound analogue of another result due to Kwerel [112] is obtained: 
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Corollary 3.7. Let A lr ...,A n be a finite collection of events where n>3. Then, 
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The following corollary generalizes Corollary [33 Its upper bound analogue [[3j[4]] is 
related to an inequality of Galambos and Xu [7j . 

Corollary 3.8. Let A 1 ,...,A„ be events. Then, for m = 0,...,n-l we have 
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Proof. For any M c {1,. . ., n) let Gm denote the join of the complete graph on M and the 
edgeless graph on the complement of M. It follows by induction on |M| that each Gm 
is chordal. By averaging the right-hand side of © over all Gm with \M\ = m we obtain 
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|M|=m 7 clique of G M 

Now follow the lines of the proof of Theorem 4.4.7 in [3] pp. 42-43], where it is shown 
that the bracketed terms in (11411 and (1 1 5D are equal. Thus, the corollary is proved. □ 



4 Open problems 

Normally when applying Bonferroni inequalities, joint probabilities of up to k events 
are given, and reasonable estimates on the probability of their union are wanted. To 
obtain a minimum upper bound, respectively maximum lower bound on this probabil- 
ity based on the inequalities in ([2]-[5]), a chordal graph having clique number k or less is 
sought which minimizes, respectively maximizes the bound on the right-hand side of 
these inequalities. This seems to be an intrinsic computational problem. Even in the 
particular case of a tree, where Hunter's method [9] gives a polynomial time algorithm 
for obtaining a minimum upper bound, no efficient algorithm is known for maximizing 
the bound in (II It . We leave these computational problems for future research. 
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Appendix: Application to Network Reliability 



Consider the directed network N in Figure 2(a)| in which all nodes are perfectly re- 



liable, and the arcs are subject to random and independent failure where all failure 
probabilities are assumed to be known. We are interested in lower bounds on its source- 
to-terminal reliability R st (N), which is the probability that there is a directed path from 
s and t consisting of non-failing arcs only. Obviously this is the case if all arcs in the di- 
rected paths 15, 136, 245, 26 are operating. For this ordering of s, f-paths, let A,- denote 
the event that all arcs in the ith path are operating (i = 1,...,4). Then, we have 

R st (N) = Pr(A 1 UA 2 UA 3 UA 4 ). 
If we assume that the reliability of each arc is p, then by the sieve formula, 

Rst(N) = 2p 2 + 2p 3 - 5p 4 + 2p 5 . (16) 

Corollary [3J] (for G = CD — CD — CD — ©) and Corollary I3.4l give 

R-st(N) > p 2 + p 3 - p 4 -ip 6 , (17) 
RstiN) >p 2 + p 3 -fp 4 - |p 6 . (18) 

In order to compare these bounds with the classical Bonferroni bound of degree 2, we 
put r = 1 in (TJJ. In this way, we obtain 

R st (N) > 2p 2 + 2p 3 -5p 4 -p 6 . (19) 



Figure 2(b) shows the reliability function ()16[) and the lower bounds ([17H19D for any p 
in the unit interval [0,1]. It turns out that for large values of p, the bounds in (1171) and 
()18[) are closer to the exact reliability value than the classical Bonferroni bound. 

From a practical point of view, larger values of p are more relevant than lower values 
of p, since network components are usually highly reliable. 




(a) Network N with terminal nodes s and t. (b) Lower bounds for R st (N). 

Figure 2: Application to network reliability analysis. 
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